A LIMIT THEOREM FOR PARTICLE CURRENT IN THE 
SYMMETRIC EXCLUSION PROCESS 



ALEXANDER VANDENBERG-RODES 

Abstract. Using the recently discovered strong negative dependence proper- 
ties of the symmetric exclusion process, we derive general conditions for when 
the normalized current of particles between regions converges to the Gaussian 
distribution. The main novelty is that the results do not assume any transla- 
tion invariance, and hold for most initial configurations. 



1. Introduction 

The exclusion process on a countable set S is a continuous-time Markov process 
describing the motion of a family of Markov chains on S, subject to the condition 
that each site can contain only one particle at a time. With the assumption that 
the jump rates from sites x to sites y satisfy p(x, y) = p(y, x), the resulting process 
is termed the symmetric exclusion process (SEP). See |14j for the construction and 
the general ergodic theory. 

In conservative particle systems such as the exclusion process and the zero-range 
process - systems where particles are neither created nor destroyed - one topic of 
study is the bulk flow or current of particles. By this we mean the net amount of 
particles that have flowed from one part of the system into the other. Finding the 
expected current in such systems is usually quite straightforward, however, given 
the interdependence of the particle motions, characterizing the current fluctuations 
is a harder problem. In the case of asymmetric exclusion on the integer lattice, the 
variance of the current as seen by a moving observer has been shown to have the 
curious order of £ 2 / 3 , with connections to random matrix theory [3ll8l[21]. 

For symmetric exclusion on Z, when only nearest-neighbor jumps are allowed, 
the current flow is intimately tied to the classical problem of determining the mo- 
tion of a tagged particle. This is especially clear when the process is started from 
the equilibrium measure v p - the homogeneous product measure on {0, 1} Z with 
density p. Since particles cannot jump over each other, and the spacings between 
subsequent particles are independent geometric- (p) random variables, the tagged 
particle's displacement is asymptotically proportional to the current across the ori- 
gin. In this case Arratia [2] gave the first central limit theorem for a tagged particle, 
and Peligrad and Sethuraman [18] showed process-level convergence of the current 
(and hence of a tagged particle) to a fractional Brownian motion. Non-equilibrium 
results were obtained by Jara and Landim [11U12] under a hydrodynamic rescaling 
of the process, even with non-translation invariant jump rates (quenched random 
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bond disorder) . The heat equation machinery used there requires the initial distri- 
butions to be smooth profiles, giving results only in an average sense. 

More recently, Derrida and Gerschenfeld [BJ applied techniques used for the more 
difficult asymmetric exclusion [21) to SEP, obtaining the asymptotic distribution 
of the (non-normalized) current. Although their results are sharp, translation in- 
variance of the jump rates and a step-initial condition seem to be required by that 
approach. 

Meanwhile, a general negative dependence theory with application to the sym- 
metric exclusion process was developed by Pemantle [19] , and Borcea, Branden and 
Liggett [4], which had immediate application to the current when SEP is started 
from a deterministic initial state [16] . 

In this paper wc further exploit the negative dependence theory in this direc- 
tion, obtaining a central limit theorem for the current throughout a wide range of 
transition rates and initial conditions. 

2. Particle current 

The original problem as described by Pemantle |19j was proved and generalized 
to the following by Liggett [IB] . Consider SEP on Z with translation invariant 
transition probabilities that describe a random increment with finite variance, i.e., 

n 2 p(0, n) < oo. 

n>0 

Start with particles initially occupying the whole half lattice {i£Z;i<0}. Then 
the current of particles across the origin after time t, 

W t = J2vt(x), 

x>0 

satisfies the central limit theorem 
W t - EW t 

— ; 7V(0, 1) in distribution. 

It was conjectured in [TBJ that this result would also hold in the case where the 
transition probabilities lie in the domain of a stable law of index a > 1. We will 
show this in Section 5. 

In this paper we consider the following general setting: Let S be an arbitrary 
countable set. For a partition 5* = A U B, we think of the net current of particles 
from A to B to be 

w t = w+(t)-w-(t), 

where W + {t) is the number of particles that start in A and end up in B at time t, 
and W~ (t) is the number of particles that start in B and end up in A. As the usual 
construction of SEP does not distinguish particles, we make this quantity rigorously 
defined through Harris' "stirring" representation, as used by De Masi and Ferrari 
[5]. The key to the stirring representation is to notice that in SEP, particles and 
holes both have the same transition rates. Hence we first define a larger process of 
randomly jumping labels, which are later reduced to either a particle (1), or hole 
(0). 

At time t = 0, we place at each site ieS the label x. For each unordered pair 
{x,y) of points, we place a Poisson process (clock) N x > v with parameter p(x,y). 
When the clock N x,v rings, the labels at x and y switch places. Let £f denote 
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the position at time t of the label x, so in particular £q = x. Under reasonable 
conditions on the rates, the random process {£f;£ G S, t > 0} is well defined on a 
set of full measure. Let L t (x) denote the label occupying site x at time t. Given 
an initial condition r\ G {0, we set 

Vt(x) = r)(L t (x)). 

Notice that if the clock N x ' y rings at time t, this produces an effect on the state rjt 
if and only if there is one particle and one hole between sites x and y\ in that case 
they switch locations. This gives one construction of SEP. 
Define the current from A to B as 

(!) W t = E vi^htfeB} - E vWhtfeA}- 

x£A x£B 

This is well defined for any initial condition rj as long as 

(2) 4H l ^ t eA } ) <oo. 

When r\ contains only finitely many particles, the current can be written as just 

Wt = E {^t( x ) ~^ x )}- 

xeB 

This coincides with the definition given above, because 

E \V( L t( x )) ~ »?(»)} = EE J 7(y) 1 {£ t w =a:} - E 

i£B xeB y£S yeB 

= E ^^eii} - E ^ 
yes a es 

= E ^(^{tfefi} - E ~ 1 {€?es}]> 

a eA y£B 

which is precisely the expression (fTj). 

Recalling that •) give the (symmetric) transition rates of individual particles 
in the exclusion process under consideration, we henceforth let X t be the one- 
particle Markov chain on S with those transition rates. 

For instance, suppose S = Z and A — {x < 0}, B = {x > 0}. Each £f has the 
same distribution as X t started from the site x, though of course for different x 
the Markov chains are highly dependent. Then under very mild conditions on the 
rates, such as 

p(x,y) <C\x-y\- a ~\ a>l, 
we can compare (using a coupling argument) X t to a translation-invariant random 
walk Zt, having finite first moment, to show that P x (X t < 0) < P°(Z t > x). 
Condition © then holds, because 

E (E !{€f <o}) < E P °( Z * ^ x ) = ^ Z t) < 

^a:>0 ' a:>0 

We say that the partition S — AUB is balanced if there is a c > 0, not depending 
on x, such that 

(3) c < liminf P x (X t E A) < lim sup (X t G A) < 1 - c. 

* — > ' 00 t— s-oo 

Here is our main theorem: 
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Theorem 1. Let S = A U B be any balanced partition of S , and rj 6 {0, 1} S be 

a (deterministic) initial condition for r\t - the symmetric exclusion process on S. 
Suppose JJ|) holds at all times, and that 

(4) supE"f V {l-rh(x))\ =oo. 

V r,(,) = l ) 

Then the current of particles between A and B satisfies the central limit theorem 

W? - EWT 

^/Var ny 

Furthermore, we have the following rate of convergence in the Levy metric: 
d(W?,N) < C(VarW?)-*. 

Condition (j4]) is a measure of how rigid the system is: by varying the time 
parameter, the expected number of initially occupied sites that are then empty 
needs to be unbounded. 

The reader can skip to the last section to see these conditions checked for a 
couple of examples. 

3. Negative dependence and SEP 

Because of the hard-core repulsion of particles, the Symmetric Exclusion process 
tends to spread out more than independent particles would. One example of this 
is the following correlation inequality of Andjel pQ: for disjoint subsets A,B of S, 
and starting configuration r\, 

(5) P n (ri t ee 1 on A U B) < P''(ry t = 1 on A)P n (r) t = 1 on B). 

There is already a well-developed theory of positive correlations, with results such 
as the celebrated FKG inequality. There, one states that a measure fi is positively 
associated if for all monotone increasing functions f,g — assuming the natural partial 
ordering on {0, 

fgdfi > J fdfx J gdfi. 

Many processes with spin-flip dynamics - such as the Ising and Voter models - 
are known to preserve positive association. That is, assuming an initial distribu- 
tion that is positively associated, the distribution of the process at later times is 
still positively associated. One may consider the following analogue for negative 
correlations: we say that /i is negatively associated if 

(6) J fgdfi < J fdfi J gdfi, 

for all increasing functions /, g that depend on disjoint sets of coordinates. The 
latter condition is a reflection of the fact that any random variable is positively 
correlated with itself. Unfortunately SEP does not preserve negative correlations 
[15] . however, there is a useful subclass of measures - introduced in [4] - that is 
preserved by the evolution of SEP. 

A multivariate polynomial / G C[zi, . . . , z n ] is called stable if 

(7) Jm(zj') > for all 1 < j < n => f(zi, z n ) ^ 0. 
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We then say that the probability measure fi on {0, 1}" is strongly Rayleigh if its 
associated generating polynomial 

is stable. 

In the setting {0, 1} for 5 infinite, we say that the measure \x is strongly Rayleigh 
if every projection of \i onto finitely many coordinates is strongly Rayleigh. It is 
easy to check that product measures on {0, l} 5 are strongly Rayleigh. 

Two key results (among many) were shown in [4]: 

(1) Strongly Rayleigh measures are negatively associated. 

(2) The evolution of SEP preserves the class of strongly Rayleigh measures. 

The distributional limits found in [16] relied upon the following result: 

Proposition 1. Suppose /i is strongly Rayleigh and T C S. Then X^gr 7 ?^) 
has the same distribution as ^2 x€ tCx, for a collection {Cx',x G T} of independent 
Bernoulli random variables. 

Combining the above proposition with standard conditions for convergence to 
the normal distribution yields a central limit theorem for strongly Rayleigh random 
variables. 

Proposition 2. Suppose that for each n the collection of Bernoulli random vari- 
ables {r] n (x);x £ S} determines a strongly Rayleigh probability measure. Further- 
more, assume the variances Var(^ 2 , gS rj n {x)) — ¥ oo as n — >■ oo. Then 



JV(0,1) 



as n 



These kinds of results were already known for determinantal processes - see 
[9lfT7] - although this is hardly a coincidence, as a large subset of determinantal 
measures are strongly Rayleigh [U Proposition 3.5]. 

During the writing of this paper, it became apparent that the argument given 
in [16] for the first proposition only holds for finite subsets T. While the results 
here and in [16j can be modified to accommodate this deficiency, in section 4 we 
will give a proof for the infinite case that may be of independent interest. 

The proof of Theorem [T] hinges upon the following proposition: 

Proposition 3. Under the conditions of Theorem^ 

Var Wt — > oo as t — > oo . 

Suppose we start with r\ = 1 on A and r\ = on B. It is trivial to check that 
product measures on {0, 1} S are strongly Rayleigh, so by Proposition [5] in order to 
show convergence to the normal distribution we need only show that the variance 
Var Wf — > oo as t — > oo. Then we may write the current variance as 

Var(W?) = Var ( ^ r, t (x) \ - £ Var(^(x)) + £ Go^ t {x)Mv))- 

^iEB ' x£B x,y£B 

x=£y 

Because of the negative association property of symmetric exclusion, all the off- 
diagonal covariances (x ^ y) are negative, while of course the diagonal terms are 
positive. The approach in [16^ for the Pemantle problem described above was to 
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compute the exact asymptotics of the diagonal terms, and then estimate the neg- 
ative (off-diagonal) terms to be at most some fixed percentage smaller. Getting 
tight-enough bounds on the negative terms was already tricky in that case - con- 
sidering even slightly more general transition functions p(x, y) seems to require 
quite delicate analysis to obtain bounds that even approached the positive terms' 
asymptotics. To get around this obstacle, we do a generator computation in or- 
der to rewrite the positive variances, and obtain term-by-term domination of the 
off-diagonal covariances. 

There is one additional complication when extending the result to the more 
general initial conditions of Theorem [TJ when n contains infinitely many particles 
in both A and B we cannot write as a convergent sum of occupation variables. 
Instead, we approximate by considering initial conditions with only finitely many 
particles. Consider first the following lemma. 

Lemma 1. Suppose S n /*• S is an increasing sequence of finite subsets, and for 
r\ £ {0, 1} S define r/ n (x) = Then for fixed t > such that fl|) holds, 

W? -> Wt in L 2 as n — ^ 00. 

Proof. Using the stirring representation ([T]) and the inequality (a — b) 2 < a 2 + b 2 
for a, b > 0, 

E{W n _ W ff < E ( hSBBMx) - V n {x)jf 

x£A 

xeB 

Both expectations are dealt with identically, so we consider here only the first one. 
Expanding it gives 

£ El { £. 6J , } l {6?6lJ} (» 7 ( a - rTix))^) - rj»(y)) 

x,y£A 

^ ( E E 

(r,(x)-ri n (x))) + (n(x)-n n (x)). 
xeA xeA 

The inequality here follows by negative dependence (0, and the expectations ap- 
pearing in the last line converge to zero by Dominated Convergence. □ 

Proof of Theorem 1. For (n,r) n ) as above we note by the triangle inequality that 

(8) d{Wf, AO < d{Wf, Wf) + d(W, AO • 

Now recall from Proposition [1] that 

wf + Y,^ x )= E c& 

x£A XESn 

where the £ are Bernoulli and independent in x for each n and t. Normalize both 
sides to see that 

W= E 

x£S„ 
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Hence by Esseen's inequality [20l V, Theorem 3], 

E Var (C 



d(W? ,JV)<C 



LxES n 



E E ic - 



< C 



E Var (0 



because C Bernoulli implies that E|£ — E£| 3 < Var(£) by an easy calculation. Taking 
n — > oo above and in ([5J , and using Lemma [IJ 



d(W?,M) < C[Var(W t ")]- 
Applying Proposition [3] finishes the proof. 



□ 



Our proof of Proposition [3] relies upon the following representation for the on- 
diagonal covariances. 

Lemma 2. Let X t be defined as above Theorem^ Then for any r\ e {0, 1} S with 
finite support (i.e., rj(x) = 1 for only finitely many x), 



(9) 



EVar^z)) = f E P^,ym y v(X s )-^ x v(Xs)} 2 ds. 

lES J° x+y 



x^y 
x,y£S 



Proof. (Essentially a generator computation). From the duality theory of SEP 
[Ml VIII, Theorem 1.1], E^O) = E x r)(X s ). So we can write 

E Var %(*) = E { E V(z) - (ev^) 2 } - E { E ^(^) - I E ^(V S )] 2 }. 



xES 



x£S 



Let [7 and {U(t);t > 0} be the generator and semi-group for Xt. Changing into 
the language of generators, we have 

EVarM*)) = E {U(s)v(x) ~ [U(s)v(x)} 2 }. 
Now take the derivative w.r.t. s. For the second equality below, recall that 

uf(x)= E p(*.v)[/(») -/(*)], 

for bounded /. 

d E {^'M*) - [U(s) V (x)} 2 } = E C/[C/(s)r/](x)[l - 2U(s)r)(x)} 



ds 



xES 



x£S 



E p{x,yWv{Xs)-^ x fl{Xs)][l-'M x r 1 {X s )] 

x^y 
x,y<£S 



(10) 



-J2p(x,yWv(x s ) - e x v(Xs)]{i - &v(x s ) - E x v(x s )}, 

x=£y 
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where all sums converge absolutely because 77 has finite support. Since p(x,y) — 
p(y, x), exchanging x and y in the latter sum in (|10p shows it to be its own negative, 
hence zero. We thus conclude that 

^£VarMa;))= ]T p(x, yW v (X s ) - E^(X S )] 2 . 

x^S Xj£y 

x,y£S 

Integrating from to t finishes the proof. □ 

Proof of Proposition^ Considering only initial conditions 77 £ {0 7 1} S containing 
finitely many particles, the net current from A to B can be written as 

w? = J2 - = E fa) - ^)}- 

xeB x£A 

We symmetrize this expression: 



2W? = Y,[ H ( x )Vt{%) ~ H(x)t](x)}, where H(x) 
and consider the variance: 



1 if x e B 
-1 if x e A, 



4Var(iy t I 'H^Var(7 7t (x))+ £ H (x)H(y) Cov( m (x), Vt (y)). 

x£S x^y 
x.yeS 

We first deal with the covariances above, proceeding almost identically to j!6j . 
Let {U2(t);t > 0} be the semigroup for two identical, independent Markov chains 
with symmetric kernel p{x 1 y), and let U 2 be its infinitesimal generator. Specifically, 

(11) U 2 f(x, [p(x, z)[f(z, y) - f(x, y)} + p(y, z)[f{x, z) - f(x, y)}}. 

Let {V 2 (t); t > 0} and V 2 be the semigroup and generator for that process with the 
exclusion interaction, i.e. 

(12) V 2 f(x,y) = ^2p{x,z)[f(z,y) - f{x,y)] + Y,v{y,z)[f{x,z) - f(x,y)}. 

z^y z^x 

By a slight abuse of notation, define 

T](x, y) = r)(x)r)(y) and H(x, y) = H(x)H(y). 
By duality and the integration by parts formula, 

-J2 H &V) Cov( Vt (x),ri t (y)) = *EH(x,y)[U 2 (t) -V 2 (t)} V (x,y) 

x^y x^y 

(13) =/ Y^H{x,y)V 2 {t-s)\U 2 -V 2 }U 2 {s)r } {x 1 y)ds. 

J ° **y 

Now from fTTj) and (|12|) we have that 

[U 2 - V 2 ]U 2 (s)r]{x, y) = p(x, y){u 2 {s)r)(x, x) + U 2 (s)r](y, y) - 2U 2 (s)r]{x, y)} 
= p(x,2;)[E^(X s )-E^(X s )] 2 , 
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which we substitute into (fT3|). also using the fact that Vzft — s) is a symmetric 
linear operator on the space of functions on {(x, y) € S 2 ; x ^ y}: 

-J2H(x)H(y)Cov(ri t (x), m (y)) 

x=£y 

= j y^p{x,y)[&r } (X s )-W c i 1 {X s )] 2 V 2 {t-s)H{x,y)ds 
Jo , 

x^y 

< f -E x V (X s )] 2 U 2 (t - s)H(x,y)ds, 

by a standard inequality comparing interacting and non-interacting particles |141 
VIII, Proposition 1.7]. 

Combining Lemma [2] and the above estimate we obtain: 

(14) 4Var(W?)> / j/)[E^(X s ) - E x v (X s )} 2 q ^ s (x,y)ds, 

Jo / 

x=iy 

where 

q s {x,y) = 1 - U 2 (s)H(x,y) = 1 - [1 - 2P X (X S E A)][l - 2P y (X s £ A)}. 

Notice that there is a constant c' > 0, depending only on the c in (J3J, such that 
Qs(x, y) > c' for each x,y £ S and then s large enough. So for fixed T > 0, applying 
Fatou's Lemma twice, then using Lemma [5] again, 

41iminf Vax(W?) > [ V p(x, y^E^iX,) - E x r)(X s )} 2 liminf 

t— >oo — ' £— »oo 

^>(x,y)[E^(X s ) -E^(X s )]Vds 



q t - s (x,y)ds 

x=jty 

T 

I . 

c! 2_, Var(7yT(a;))- 

xes 



Now for any finite S' C S, 

Var(?7 T (a;)) = ^ ^ v(v)Pt{x, y) I -^r]{z)p T {x,z) 
xes' xeS' yes L zgs 

= Yl PT(y,x)p T {x,z) -> ^ ^ P2T{z,y), 

as S" 5, by monotone convergence. By duality again, this is precisely 

E"(l-r? 2T (y)). 

Hence by (jH), 

liminf Var(W t ' J ) = oo, 
as desired. □ 
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4. Sums of strongly Rayleigh random variables 



In this section we prove Proposition 1 for infinite-dimensional ji. By a general- 
ization of the Borel-Cantelli lemmas [7J, applied to the negatively dependent events 

r](x) — oo, /i-a.s. if E M T)(x) = oo. 
cceT ieT 
In this case the proposition is trivially true, hence we may assume that 

r)(x) < oo. 



First take an increasing sequence of finite subsets T ri 
the polynomials 

Q n (z) =Wz^T» r]{x) . 

The limit 

Q(z)=Wz^tv(x) 
exists, and in fact Q n — > Q uniformly on compact sets. Indeed, 



T. Now define for z G 



\Q n (z)-Q(z)\<W 



Tf(x) 



1 — z 52xer\r, 



< E^maxjl, |z|E- € t "^}l { ^ on t\t„} 



Now for \z\ = r > 1, r is an increasing function of t, so the negative dependence 
property of \i ([6]), implies that 

(15) E^r^^T •»(*) < Y[ EV (s) = JJ [1 + (r - l)W-q(x)] < e r ^ < oo. 

(The last inequality uses the estimate 1 + x < e x .) Dominated convergence then 
gives the normal convergence Q n —¥ Q. In particular, Q is entire. 

By stability ([7]). Q n (z) has only real zeros, furthermore, the zeros must all be 
negative because the coefficients of Q n are all non-negative. By classical theorems 
on entire functions 13, VIII, Theorem 1], the limit Q(z) has the form 



Q(z) = Ce"" JJ 



k=l 



1 



2 



for some a < 0, and < with X) l a fcl 1 < 00 ■ It is enough to show that a = 0, 
because, as Q(l) = 1 we can solve for C to obtain 

oo oo 

QW = II^ri = IIb^ + (i-p*)] ) 

fc=l k k=l 

where we set pk = 1/(1— a/-). But this last expression is just the generating function 
for the sum of independent Bernoulli r.v.'s having the parameters pk- 

To obtain a = 0, we show that |Q(z)| < e c ' 2 for any c > and \z\ large enough. 
It is clear that |Q(z)| < Q(|z|), hence we consider only z = r > 1. Recall from (fl5|) 
that 

Q(r) < IJ[l + (r-l)E"r7(a;)]. 
xeT 

Let 

a, =E"t ? (o;). 
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With a — J2x<et ax < °°i n °t e that #{x : a x > r 1 / 2 } < ar 1 / 2 . Then by a trivial 
bound we have 

Q(r)<( n A( n e( "^) 

< r^ 2 exp ((r - 1) £ a A 
As r —> oo the sum inside the exponential goes to zero, which concludes the proof. 



5. Examples 

1. Consider S = Z, partitioned into A = {x < 0} and £? = {x > 0}, with 
translation invariant rates p(0,x) in the domain of a symmetric stable law with 
index a > 1. That is, 

y>x 

for a slowly varying function L. Consider the step initial condition rj with particles 
at all x < 0. The balance condition holds by the central limit theorem for random 
variables in the domain of attraction of a stable law. By duality and translation 
invariance, 

x<0 x<0y>0 n>0 

(In fact, it grows at rate t 1 /"). This shows ([3]), and the above expression is the 
same as WWt, so all conditions of Theorem [T] have been verified. 

2. Now consider the one-dimensional exclusion process in an random environ- 
ment, with the same partition as in the previous example. The random environment 
is described by {u>i}, an iid family of random variables with uj € (0, 1] almost surely 
and E— < oo. For each realization (. . . , uj-\, (da, u)\, . . .), we consider the exclu- 
sion process with the rates p(i, i + 1) = p(i + 1, i) = Wj. By the result of Kawazu 
and Kesten [10], we know that the process {X n 2 t /n} converges weakly to a scaled 
Brownian motion, from which follows the balance condition. By the remarks after 
@ we know that the current has finite expectation for any initial placement of par- 
ticles. Let us consider the case where we pick a realization r\ of the homogeneous 
product measure v p . Then almost surely- v p there are infinitely many x £ Z such 
that r\{x) — 1 and r\{x + 1) = 0, so from the ergodicity of the environment, 

co= Y. P x {X t =x + l) < P x (r,(Xt)=0), 

rj(x) — l rj(x) — l 

>)O+l)=0 

which shows by duality. 

Acknowledgments. This article is part of the authors' thesis under T. M. Liggett, 
whom the author would like to thank for his advice and encouragement. 
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